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Abstract 

We investigate numerically by a conservative difference scheme in complex arithmetic the head-on 
and takeover collision dynamics of the solitary waves as solutions of linearly Coupled Nonlinear 
Schrodinger Equations for various initial phases. The initial conditions are superposition of two 
one-soliton solutions with general polarization. The quasi-particle behavior of propagating and 
interacting solutions in conditions of rotational polarization is examined. We find that the total 
mass, pseudomomentum and energy are conserved while the local masses, individual and total 
polarization depend strongly on the linear coupling and the initial phase difference. We also find 
out that the polarization angle of the quasi-particles can change independently of the interaction. 
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1. Introduction 

The system of Coupled Nonlinear Schrodinger Equations (CNLSE) is an attractive because of 
several reasons. From among them: CNLSE are a soliton supporting system, the system is a proper 
"a spring board" for investigate and track the quasi-particle (QP) behavior of the solitary waves; it 
possesses very rich phenomenology which can be extended adding new terms in the equations; the 
vector structure of the system and their solutions admits to study and treat the polarization vector 
lfl6l . 0. The initial polarization of CNLSE and its evolution during the QP interaction is a 
very important element. There are numerous experimental observations (see flU and cited there 
literature) of the effects connected with this quantity. For the full fledged CNLSE with general ini- 
tial polarization, no analytic solution is available. It can be obtained, e.g., via adequately devised 
numerical scheme. An implicit scheme of Crank-Nicolson type was first proposed for the single 
NLSE in [8J. In order to adapt it to CNSLE, internal iterations were applied in Q. The latter 
yielded two very important advantages: fully implicit scheme and implementation of the conser- 
vation laws on difference level within the round-off error of the calculations. The above scheme 
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was extended to complex arithmetic in [10] where the computer code for Gaussian elimination 
with pivoting of [2] was generalized for complex- valued multi-diagonal band algebraic systems. 
Having in mind the importance of initial polarization, we investigated in [1 1 J the collision dynam- 
ics for circularly polarized solitons based on sec/i-functions and found out that there exist infinite 
number of Manakov two-soliton solutions preceded by polarization discontinuity (shock) on the 
place of interaction. In order to extend the range of investigations in the case of general polariza- 
tion we established an auxiliary conjugate system of nonlinear ordinary differential equations [fT2l 
in order to generate numerically initial elliptically polarized soliton solutions and applied them for 
numerous head-on [fT3l and takeover collisions 0. To enrich the concept of the polarization we 
involved into considerations linear coupling of sec/i-like solitons which is a generator of their ro- 
tational polarization (breathing) and/or gain/self-dissipation [14]. In this work we aim to conduct 
series of simulations and to track the particle-like behavior of the interacted localized non-sech 
waves with arbitrary (elliptic) initial polarization. 

2. Problem formulation 

CNLSE is system of nonlinearly coupled Schrodinger equations (also called the Gross-Pitaevskii 
or Manakov-type system): 

iV>t = PAx + NVf + ("i + 2a 2 )\<f)\ 2 }ijj, (la) 
i<t>t = P4>*x + KM 2 + («i + 2 «2)M 2 ]0, (lb) 

where (3 is the dispersion coefficient, aii describes the self-focusing of a signal for pulses in bire- 
fringent media, a 2 (called cross-modulation parameter) governs the nonlinear coupling between 
the equations. When a 2 = 0, no nonlinear coupling is present despite the fact that "cross-terms" 
proportional to oi\ appear in the equations and the both equations admit identical solution equal 
to the solution of single NLSE with nonlinearity coefficient a = 2a.\. Now we turn our consid- 
erations to a modified CNLSE where the nonlinear coupling is trivial but is replaced by so-called 
linear coupling. Nevertheless the system remains nonlinear 

tyt = M** + «i(M 2 + |0| 2 )^ - T0, (2a) 
i(pt = ^ xx + a 1 (\^\ 2 + ^\ 2 )^-T^. (2b) 

The magnitude of linear coupling is governed by the complex- valued quantity T. 3?r governs 
the oscillations between states termed as breathing solitons and actually seems their frequency 
of oscillation, while QT describes the gain/self-dissipation behavior of soliton solutions. As it is 
shown in the breathing is present free of the interaction. The relation between the Manakov- 
type system ([[]) (when a 2 = 0) and the Linearly CNLSE ([2]) is given by the substitution (for details 
see and O) 

$ = $ cos(rt) + i$ sin(rt) , = $ cos(H) + sin(rt) , (3) 

where $ and \I> are soliton solutions of ([T]). These properties and references ([3]) of the solutions 
allow and argue us to solve (|2]) instead ^ in order to get more information about the breathing 
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behavior of the soliton solutions. Let us remind that in ffT4l we considered the problem in question 
but in the particular case when <3> and \I/ were assumed to be certainly sec/z-solutions of ([T]). 

In this paper the solutions are pulses whose modulation amplitude is of general form (non- 
sech) and their polarization^] rotates with time. This determines the choice of initial conditions 
for numerical investigation of temporal evolution of interacting solitons. As usual we concern 
ourselves with the soliton solutions which are localized envelopes on a propagating carrier wave. 
We assume that for each of the functions 0, ip the initial condition has the general type 

x(x, t; X, c, n x ) = A x (x + X — ct) exp {i \n x t — \c{x — X — ct) + <5 X ] } (4) 

where x stands for (tp,(f>), c is the phase speed of the envelope, X is the initial position of the 
center of the soliton; n(n^, n^) is the vector of carrier frequencies of the components; 5(5^, 5$) 
is phase vector of the two components. Note that the phase speed c must be the same for the 
two components ip and <\>. Otherwise they will propagate with different phase speeds and after 
some time the two components will be in two different positions in space, and will no longer form 
a single structure. For the envelopes (A^,A^), 9 = arctan(max |0|/max \ip\) is a polarization 
angle. We are interested in solutions with initially nontrivial carrier frequencies (n^,, n^) ^ (0, 0), 
i.e., with elliptic or circular polarizations. Let us emphasize that these qualities are intrinsic to non- 
sech-\\ke solutions. Having in mind ([4]) we implement straightforward manipulations yielding 

[A'lp + (rty + \c 2 )A^ + a x (A\ + Al) A J exp {i [n^t - \c(x - X - ct) + S^\ } cos(rt) 

+i [A'j, + (t^ + \c 2 )A<p + ai (Al + Al) Af] exp {i [n^t - \c(x - X - ct) + fy] } sin(rt) = 0, 

(5a) 

[A'i + (n* + \c 2 ) A^ + a x (A\ + A\) A J exp {i [ n<j> t - \c(x - X - ct) + 8A } cos(rt) 

+i [Al + (n f + \c 2 )A^ + a x (A\ + A\) A^\ exp {i [n^t - \c(x - X - ct) + fy] } sin(rt) = 0. 

(5b) 

The above two equations are true if and only if the two components of the envelope are governed 
by the following coupled system of differential equations: 

Al + {jH + \c 2 )A^ + ai (A\ + A\) A^ = 0, (6a) 

A 'i + K + A <P + «i ( A l + A l) A <t> = °- ( 6b ) 

We find out that this system is explicitly equivalent to the bifurcation conjugate system^] corre- 
sponding to Manakov case with a 2 = (see [|T3l and lfT2ll ). Here we aim to trace the strong linear 
coupling on the interaction in conditions of general initial polarization. After that we construct the 
initial conditions for x from as a superposition of two two-component solitary waves located 
at positions Xi and X r and propagating with phase speeds q and c r 

x(x, 0) = xi(x, 0, c h n h X h Si) + x r (x, 0, c r , n r , X r , S r ). (7) 

The initial distance \Xi — X r \ is supposed to be large enough in order the solutions not to overlay 



'The latter determines one more kind of polarization called rotational polarization. 

2 The system admits bifurcation solutions since the trivial solution obviously is always present. 
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(c) {5} = 180° 

Figure 1: Head-on collision of QPs with initial circular polarization for n^, = n r ^, = ni^ = n r $ = — 1.5, q = 
— c r = 1, a>i = 0.75, T = 0.175, Xi = —40, X r = 40 and different initial phase differences. 

in the time moment t = 0. The system ([6]) is solved numerically and their solutions are readily 
used as an approximate initial condition for the unsteady computations. Our aim is to understand 
better the influence of the initial polarization and the initial phase difference on the particle-like 
behavior of the localized waves (QPs). They survive the collision with other QPs (or some other 
kind of interactions) without losing its identity. The initial difference in phases can have a profound 
influence on the polarizations of the solitons after the interaction as well as on the magnitude of 
the full energy and the amplitude of breathing. 

Before turning to the numerical investigation we mention here that the system Eq. ([2]) possesses 
three conservation laws when asymptotic boundary conditions are imposed, namely when ■0,0—)- 
for x — > ±oo. Following (3113 we define "mass", M, (pseudo)momentum, P, and energy, E as 
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Figure 2: Mass dynamics corresponding to Figure[T] 
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Figure 3: Influence of the initial phase difference on the energy dynamics corresponding to Figure [T] 
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I W + M 2 )d£, P = - / + H x )dx, E = I Hdx 
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Li 



+ - 2r[^)] 



(8) 



is the Hamiltonian density of the system. Here —L 1 and L 2 are the left end and the right end of 
the interval under consideration. The conservation laws read 



dM 



0. 



dP 
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0. 



dE 
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Figure 4: Individual and total polarization dynamics in dependence on initial phase difference corresponding to 
Figure [T] 



which means that for asymptotic boundary conditions CNLSE admit at most 3 conservation laws, 
i.e., the system ([2]) is non-fully integrable. 

3. Numerical method 

In order to be able to obtain reliable results for the time evolution of the solution, one needs to 
devise a difference scheme that represent faithfully the above conservation laws. Such a scheme 
was proposed in 0, and applied in 0. This scheme was based on a fast Gaussian elimination 
solver for multi-diagonal systems Q. Consequently, the above mentioned scheme was imple- 
mented for complex arithmetic in iTTOll . where a complex arithmetic algorithms was developed 
to generalize the one from [2]. The complex-arithmetic algorithm is four times faster, and we 
will use it also in the present paper. Thus, for solving Eqs. ([2]) with the initial conditions (Q, ([7]) 
numerically, we use an implicit conservative scheme in complex arithmetic: 

? /; n + 1 _ J,n ft 



At 2h 2 



+ + (IC + T + ICI 2 ) + (IC + T + ICI 2 ) -o(C +1 + C)> dOa) 
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kn+l 
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+ ICI 2 ) + (IC + T + ICI 2 )J - g (^ n+1 + O » ( 10b ) 



on mesh (xj,r n ) with = —L\ + ih, h = (L 2 + Li)/m, i = l,...,m and r n = nAr, 
n = 0, 1, 2... It is not only convergent (consistent and stable), but also conserves mass, pseudomo- 
mentum, and energy, i.e., there exist discrete analogs M n , P n and E n , for ([8]), which arise from 
the scheme (for details see [131171 [TOlO. 
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for n > 0. These values are kept constant by the scheme during the time stepping. The above 
scheme is of Crank-Nicolson type for the linear terms and we employ internal iteration to achieve 
implicit approximation of the nonlinear terms, i.e., we use its linearized implementation Q. In 
this way the order of approximation of Eqs. ( |T0| ) is 0(At 2 + h 2 ). 

The above nonlinear scheme is implemented via internal iterations as follows: 
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The iterations are repeated (stepping up the index 'A;') until convergence is reached. If the 
internal iterations are convergent, the scheme in full steps Eqs. ( fit)] ) is absolutely stable due to its 
conservative nature. For not very large time steps the iteration requires less than half a dozen loops 
to converge. This is a small price to pay having in mind the inextricably coupled five-diagonal 
complex structure of the matrix. 

The time-stepping scheme needs initial conditions of the type of Eq. Q. In the general case 
of elliptic initial polarization we discretize the system Eqs. (|6]) in the same manner as the original 
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Figure 5: Head-on collision of QPs with initial elliptic polarization for 9 = 23°44', ni^ = n r ^, = —1.1, ni^ = n r 
— 1.5, c; = — c r = l,ai= 0.75, r = 0.175, X/ = —40, X r = 40 and different initial phase differences. 



evolutionary system, and use Newton's method and hermitian splines to get approximately the 
solution [12J. The result are used as initial conditions for the problem under consideration. 

The above presented scheme and algorithm have been verified for different girds and time 
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Figure 6: Mass and pseudomomentum dynamics for [J] = 0° corresponding to Figure|5] 
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Figure 7: Influence of the initial phase difference on the total energy dynamics corresponding to Figure|5] 



increments and the approximation has been confirmed. 



4. Results and discussion 

4.1. Circular polarization. Head-on collision 

This is a special elliptic polarization with polarization angle 9 = 45°. The initial configuration 
is generated from the auxiliary bifurcation system ([6]). Because the parametric space of the prob- 
lem is too big to be explored in full, we choose = n r ^ = = = —1.5, q = — c r = 1, 
«i = 0.75, T = 0.175 and focus on the effects of [5] = 8 r — 8\. We observe that the breathing 
does not interfere with the soliton collision since they are distinguished from previously obtained 
breathers lfi~5l . HI. One sees the "breathing" of the pulses even without any interaction. This is the 
manifestation of the rotation of the polarization. Thus Y is responsible for the exchange of wave 
mass between the modes and we call it "cross" dispersion of the signals (see Figure [2]). Recall 
that these solutions apply specifically to the case a 2 = 0. Our initial condition is valid only for 
this case. In Figure [T] we present these features for fixed set of parameters varying only the initial 
phase difference {5} . As we discussed in 031 there is no velocity shift on the place of interaction. 
The individual masses start to oscillate since the same beginning with an equal period and with 
a phase shift and an amplitude depending on the concrete choice of {5} The total mass, how- 
ever, is perfectly conserved - it is equal to 6.0683. Due to the symmetry the pseudomomentum 
P = 10~ 3 -r 10~ n and evidently approximates 0. The result is the same for all considered phase 
differences (0°, 90°, 180°) (Figure [TJ. What happens with the energy one can see reviewing the 
next Figure [3] The magnitude of the energy depends essentially on the magnitude of initial phase 
difference value [5] . The dependence is non-monotonous: E = —3.078 for [5J = 0°; E = —4.15 
for [5] = 90°; E = —1.139 for [5] = 180°. Figure [3] demonstrates cogently the conservation 
properties of the used difference scheme - one observes excellent conservation of the total energy 
of the interacting soliton system during the whole time of evolution. The next Figure [4] elucidates 
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Figure 8: Individual and total polarization dynamics in dependence on initial phase difference corresponding to Fig- 
ure [5] 



the polarization evolution. The individual polarizations oscillate depending on both the magnitude 
of linear coupling Y and initial phase difference \S\ . These oscillations do not depend on the in- 
teraction but the period depends on Y. The magnitude [<5] influences the amplitude of oscillation. 
In contrast to masses (see Figure [2]) the total polarization is not constant. It breathes also though 
with very small amplitude. It is well seeing that on the place of interaction one observes a shock 
of both the individual and total polarization. The last result can be interpret as a generalization of 
the polarization conservation law for the case of rotational polarization - it breathes in the period 
but it is conserved in the whole period. 

Our observations show that these interactions are independent of initial phase shift and no- 
phase velocity shift occurs on the place of interaction. The trajectories of QP-centers suffers slight 
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Figure 9: Takeover collision of QPs with initial elliptic polarization: polarization angles Oi = 23°45', 6 r = 25°52', 
ci = 1, c. r = 0.8, niij, = n rt f, = —1.1, = n r< p = —1.5, ot\ = 0.75, T = 0.175, Xi = —20, X r = and different 
initial phase differences. 

shift. They were investigated and discussed in 031 for sec/z-profiles and initial linear polarization 
and we do not concern them again. 

4.2. Elliptically polarized solitons. Head-on collision 

In order to extend our considerations we conduct series of experiments with general initial 
polarization and nontrivial linear coupling. For convenience we use the same magnitude Y = 
0.175. On Figure |5| the case of initial elliptic polarization is present for initial phase shifts [5] = 
0°,90°, 135°, 180°. Our conclusion is that the linear coupling combined with the initial phase 
difference [5] can shift only the soliton trajectories and keeps the phase velocities unchanged. 
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Figure 10: Mass and pseudomomentum dynamics for [(5] = 0° corresponding to Figure [9] 



A nontrivial phase shift is possible after involving nonlinear coupling a 2 ^ BH. We consider 
two solitons with equal initial elliptic polarization angles of 23°44'. The initial configuration is 
generated from the auxiliary bifurcation system ([6]). The above angle corresponds to the parametric 
set riij/j = n r ^ = —1.1, = n r ^, = —1.5, C\ = —c r — 1, a.\ — 0.75. Obviously the results are 
qualitatively the same as in the previous case of circular polarization and they do not require a 
detailed discussion. On Figure [7] is plotted the energy depending on the initial phase difference. 
The correlation between them is essential. The conservation is again excellent: [<5] = 0° - 
E = -0.262; [5] = 90° - E = -0.821; [5] = 135° -E — -0.206; {5} = 180° -E — 0.640. The 
breathing behavior of the total polarization is well visible as well as the affect of the initial phase 
difference upon the amplitude of the individual polarization angles which result to the magnitude 
of the total polarization angle. Figure [5] answers convincingly the question about the dynamics of 
the rotational polarization when the type of initial polarization is general. The net conservation law 
is generalized as net breathing (oscillating) law which conserved the total polarization within one 
time period. This result has not an analytic analogue and enriches the phenomenology of CNLSE. 
Among the three conservation laws about the mass, pseudomomentum and energy we found after 
an extensive set of numerical investigations one more net conservation law. 
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Figure [9] 



Individual and total polarization dynamics in dependence on initial phase difference corresponding to 
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Figure 12: Influence of the initial phase difference on the total energy corresponding to Figure [9] 

4.3. Elliptically polarized solitons. Takeover interaction 

The pattern in takeover interaction is the same qualitatively but the required difference in the 
velocity magnitudes bears bigger trajectory shift of the slower soliton and slower trajectory shift 
of the bigger soliton (see lfT4lO . Here we consider this kind of interaction in order to implement 
an interaction of elliptically polarized solitons with different polarization angles. Also such kind 
of evolution and interaction do not possess a symmetry. Our main goal is to check and confirm 
the conclusions about the net polarization realized in the previous subsection. We choose = 
Ti r ij> = —1-1, nitf, = n r< j, = —1.5, c\ = 1, c r = 0.8, a\ = 0.75, T = 0.175 and focus on the effects 
of 6. The above parametric set results in polarization angles 9i = 23°45', 9 r = 25°52' of the 
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initial envelope configuration. The approximate initial condition is generated from the auxiliary 
bifurcation system. The solitons start with different elliptic polarizations and are an object of 
rotational polarization as well due to nontrivial linear coupling T (Figure [9]). We find out that 
the initial phase differences of the components play an essential role on the full energy of QPs 
(Figure [12|). In particular, when {6} = 0° - E = 0.0673; {6} = 90° - E = -0.976; [5] = 180° 
- E = —0.657. Just in opposite, the pseudomomentum is also conserved and does not depend on 
the initial phase difference. It is not trivial because of the nonsymmetry of the QP configuration 
and its net values P = 4.08 -r- 4.05 (Figure 10). The individual masses, however, breathe together 
with the individual (rotational) polarizations. Their amplitude and period do not influenced from 
the initial phase difference (Figure 10) and are conserved within one full period of the breathing. 
The total mass is constant. Both the individual and total polarizations breathe and suffer a 'shock 
in polarization' when QPs enter the collision. The polarization amplitude evidently depends on 
the initial phase difference (Figure 11). The above quantities are conserved within one full period 
of the breathing. In all considered cases of evolution and interaction the polarization angle of QPs 
can change independently of the collision due to the real linear coupling T. 



5. Conclusion 

In this work we are aiming to complete the information about the dynamical properties of the 
total polarization influenced by the initial parameters and, in particular, by the initial polarization 
of general type as well as by the nontrivial linear coupling T and the initial phase difference [<5] . 
For the case of general elliptic polarization there is not available an exact initial condition concern- 
ing the shape of solitary envelopes. Solving an auxiliary bifurcation system of ordinary differential 
equations we construct numerically an approximate initial conditions for the full range of polar- 
ization angles between 0° and 90° and with non-sec/i shape. Let us denote that the particular cases 
9 = 0°; 90° (left and right linear polarization) when one has an exact sec/i-like initial condition 
were implemented in detail in lfi"4ll . We consider the two main interactions - head-on and takeover 
collisions of envelopes when only a real linear coupling T is present (the cross-modulation param- 
eter ct2 = 0). Adding an imaginary part in T leads to gain/blow-up of the solution accompanied 
by a violation of the conservation laws (see [Tl4l ) and we will not concern again this property in 
the present work. The main conclusions can be summarized as follows 

• The interaction conserves perfectly the pseudomomenta and energy. The pseudomomenta 
are not influenced by the initial phase difference [<5] while the magnitude of the energy is 
strongly dependent on it. The variation is nonmonotonous. 

• The total mass is also perfectly conserved but this is realized on the inner share and exchange 
of local masses of both components of QP periodically in the time. This periodicity is 
governed by the real part of the linear coupling T and generates breathing of the masses 
among with the breathing of the solitary envelopes. The mass breathing does not depend on 
the interaction and is present during the whole time of evolution. The amplitude of mass 
oscillation depends on the initial phase difference. 

• The individual polarizations breathe similar to individual masses but in opposite to them 
with a phase shift. This is the reason the total rotational polarization to be conserved within 
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one period. On the place of interaction the individual and total polarization angles suffer 
a discontinuity keeping the else functions smooth. In other words the total polarization 
breathes also and a bifurcation occurs during the cross-section of the interaction. And this 
is the main result of this investigation. The amplitudes of the individual polarizations are 
affected by the initial phase difference. 

• We found and generalized for the case of rotational polarization one more conservation law 
- for the net total polarization. So, the linearly CNLSE possess four conservation laws: for 
the total mass, for the pseudomomentum, for the energy, and for the magnitude of the total 
polarization. Let us emphasize that the last law has not an analytic analogue and it unearths 
more intrinsic information for the propagating and interacting solitons as quasi-particles. 
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